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ABSTRACT 



We present sensitive NIR (J, H and K) imaging observations toward four lu- 
minous massive star forming regions in the Norma Spiral Arm: G324. 201+0. 119, 
G328.307+0.432, G329.337+0.147 and G330.949-0.174. We identify three clus- 
ters of young stellar objects (YSO) based on surface density diagnostics. We 
also find that sources detected only in the H and K-bands and with colors cor- 
responding to spectral types earlier than B2, are likely YSOs. We analyze the 
spatial distribution of stars of different masses and find signatures in two clusters 
of primordial mass segregation which can't be explained as due to incompleteness 
effects. We show that dynamic interactions of cluster members with the dense 
gas from the parent core can explain the observed mass segregation, indicating 
that the gas plays an important role in the dynamics of young clusters. 



Subject headings: Hn regions — infrared: stars — stars: early-type — stars: for- 
mation — stars: pre-main sequence — stars: luminosity function, mass function 



1. Introduction 

While still deeply embedded, massive star clusters keep an imprint of the physical condi- 
tions in the natal molecular cloud. Massive stars form in molecular cores with t ypical radii of 



0.4 pc, molecular numbe r densities of 3 x 10 5 cm 3 and masses of 5 x 10 3 M Q ( jPlume et al. 



19971 : iFaundez et al.l l2004f ) . Radio wavelength observations indicate that massive stars are 



1 Departamento de Astronomi'a, Universidad de Chile, Camino del Observatorio 1515, Las Condes, San- 
tiago, Chile 

2 Universidad Nacional Autonoma de Mexico, Apt Postal 3-72 (Xangari), MX Morelia, Michoacan 58089, 
Mexico 



-2- 



formed in groups (jHo fc Haschickl Il98ll ; iGaray et al.l 119931) and tha t they are generally lo- 
cated in the center of massive and dense cores (jGaray et al. 2007 ); an indication of mass 
segregation. Thus, the study of young massive clusters is essential to fully understand the 
process of massive star formation. 

Stellar mass segregation in a cluster can be either a consequence of dynamical evolution, 
or of biased massive star formation towards the core center. Dynamical mass segregation 
arises from kinetic energy equipartition between the cluster members. Massive stars will 
move slower than low-mass stars and sink to the center of the cluster in a process that takes 
place in a time comparable to the relaxation time (t r ) of the cluster (jBinney &: Tremaine 
19871 ). On the other hand, primordial mass segregation takes place much faster, in a mas- 
sive star accretion timescale. During this process, ma ssive stars are create d in the center of 
the molecular cloud, where the gas den sity is highest (jBonnell et al.ll200ll ) and/or the mass 
accretion rate is high (IKrumhola 120061 ) . Several studies have shown evidence of mass segre- 
gation in young open clusters having ages of at least one order of magnitude smalle r than 
their relaxation time (iGouliermis et al.ll2004j ; iBragg fc Kenyonll2005l ; IChen et al.ll2007l ). This 
evidence suggests that mass segregation in those clusters is primordial. However, dynamical 
timescales are often calculated taking into account star-to-star interactions only, ignoring 
the gas. The later may have an impact in young clusters, where most of the mass is still 
in the form of gas. Recent studies have shown that the gaseous component may have an 
important eff ect in the dynamics of th e stars by gravitational drag, making dynamic times 
even shorter (jEscala et al.ll2003l . 12004 . in § 5.2). Thus, studies of mass segregation in young 
clusters should include the effects of dense gas in the cluster dynamics. On the other hand, 
Ascenso et al. (2009) have argued that sample incompleteness in observational data can 
mislead mass distribution analysis, producing mass segregation in non-segregated clusters. 
Therefore, to assess whether or not there is mass segregation in young clusters it is necessary 
to use robust mass segregation indicators. 

In this paper we study the spatial distribution of stars in four luminous massive star 
forming regions located in the Norma spiral arm, at distances between 5.4-7.3 kpc from the 
Sun. The regions were observed in the near-infrared J, H, and K bands at the VLT. Garay 
et al. (2009, hereafter Paper I) present dust mm-continuum observations towards these 
same regions. Target selection and data reduction are presented in § |5J Results, including 
identification of young stars, analysis of their spatial distribution and extinction maps are 
presented in § |3j In § [4] we study the radial distribution of cluster members of different 
masses. We also discuss the age of the clusters, the evidence of mass segregation and the 
presence of gaps in the color-magnitude diagrams. Conclusions are presented in § 
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Observations and data reduction 



2.1. Source Selection 



We observed the regions G324.201+0.119, G328.307+0.432 G329 337+0.147 and G330.949- 
0.174, selected from the CS(2— >1) survey of iBronfman et al.l (119961 ) by being among the 
brightest objects at far- infrared (FIR) wavelengths (L PIR > 6 x 10 5 L ). The four regions are 
located in the tangent of the Norma spiral arm (see Tabled]) and their CS(2^1) line profiles 
show broad line- widths (FWHM between 6.1 and 6.9 km s _1 ) and the presence of wings, 
making them good candidates for young clusters containing massive stars. 

The LSR velocities of G324.201 and G329.337 are close to the terminal velocities, which 
are the maximum allowed by pure circular motion along their lines of sight, thus, we locate 
them at the sub-central distances. These two objects are located well away from the centers 
of their parental molecular clouds as traced by CO (IBronfman et al.lll988l ). We assign the 
near kinematic dis tances based on the d istance of the molecular clouds in which they are 
deeply embedded ( IBronfman et al.l |2000| ) to the sources G328.307 and G330.949. These 
sources represent massive star forming regions in different environments, i.e., in the center 
and in the periphery of giant molecular clouds. 



2.2. Near-IR data 

The four sources were observed with ISAAC at the VLT-Antu in April 1999, using an 
1024 x 1024 pixel InSb Aladdin detector array with a scale of 0.147 arcsecond per pixel. 
The typical seeing during the observations was 0.4 arcseconds. We observed in the J, H 
and K bands centered at 1.24 1.65 and 2.16 /im respectively. The observations were done in 
dithering mode with 30 arcsecond shifts. The integration times are shown in Table [2j 

The data were corrected for bias, dark current, and pixel to pixel sensitivity variations. 
The XDIMSUN task in IRAF was used to create a sky image using the median time average 
of the dithered observations. We subtracted the sky image and corrected the field distor- 
tion.] by transforming the individual images to a central reference image using second order 
transformations in X and Y as recommended by ESO. This had a minor impact (< 1") in 
the source position at the edges of the detector. The field distortion correction, however, im- 
proved the roundness of the stellar images increasing the photometric precision. Finally, we 
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median combined the resulting images and defined the coordinate system using WCSTools^] 
and the 2MASS catalogue. 



2.3. Photometry 

We searched for star candidates using the DAOPHOT package in IRAF with a 3 sigma 
detection threshold. Then we performed aperture photometry using an aperture radii of 
0.45 arcseconds. The typical full width at half maximum (FWHM) of the detected sources 
was 2.7 pixels or 0.4 arcseconds. Sky emission was eliminated by subtracting the median 
value inside a ring with inner radius of 1.5 arcseconds and 1.5 arcseconds of width. Finally, 
a carefully visual inspection was done over the whole mosaic and in particular where there 
was extended emission in order to reject diffuse features and saturated stars. 

Calibration was performed by minimizing residuals to corresponding 2MASS detections. 
The RMS for residuals between the data sets used was less than 0.05 magnitudes in all filters. 
Finally, we combined J, H, and K photometry by merging the data. A source in different 
bands was identified as the same if the difference in their position was less than 0.3 arcseconds. 
Figure [1] shows the photometric error dependence over magnitude (upper panel) and over 
distance from the center of the mosaic in the K-band for G324.201 (lower panel). We only 
use sources with an error of less than 0.2 mag in all bands. This corresponds to 88% of the 
detected sources in G324.201 and G328.307, and to 77% in G329.337 and G330.949. Sources 
with an error of more than 0.2 mag are distributed uniformly over the mosaic. 



3. Results 



Figure [2] shows composite three-color images of the four regions observed. In G324.201 
the field of view is not fully covered in the H-band. Extended emission is clearly seen in all 
regions with the exception of G330.949 and shows a more complex structure in G324.201 and 
G328.307 than in G329.337. Since this emission is seen in the J and K-bands, we hypothesize 
that most of it corresponds to stellar emission reflected by the dust. Table [3] lists the total 
number of stars detected in each of the bands with error < 0.2 mag. The completeness of 
our data (shown in Table HJ) corresponds to the magnitude at which at least 90% of the stars 
is detected and it was estimated by the fraction of artificial stars recovered. Figure [3] shows 
the histogram of the detections in each band for G328.307. 



2 http: / /tdc-www. harvard.edu/software/wcstools/ 
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3.1. Color-color diagrams 



Color-color diagrams of the stars within the observed regions are shown in Figure HI 
Most of the stars are located in the region between the reddening vectors and they corre- 
spond to normally reddened stars. A small percentage of stars (between 2.5 and 10%) are 
distributed above the main sequence. This is typica lly seen in other regions of star forma- 
tion, (eg. M17, Westerlund 2 and NGC7538, from Ijiang et al.l booi kscenso et all 120071 : 
Balog et al.l 12004 respectively). We found that these sources are distributed uniformly over 
the field and do not have a characteristic brightness. The scatter in the color-color diagram 
was quantified by computing the mean distance to the giant branch, A, of the sources located 
above it. This scatter is not due to systemic errors in the reduction of the data and affects 
a small percentage of sources. 

We consider sources with NIR-excess those located at least one A to the right of the 
reddening vector (marked as asterisks in Fig. H]). Some of these sources may be background 
galaxies. To estim ate the number of background contaminants we used the K-band galaxy 
number count from lBlanc et al.l (120081 ) and found a maximum of 10 galaxies in an area of the 
same size as the area covered by our data (with K < 18 mag). Since the number of galaxies 
in the FOV is small, we assume that all the NIR-excess sources correspond to young stellar 
objects. We found 236, 120, 186 and 228 sources with NIR-excess in G324.201, G328.307, 
G329.337 and G330.949, respectively. 



3.1.1. Extinction law 



We used the NIR color-color diagrams to determine the extinction law Ej_h/ Eh-k in 
each region. In order to do this, we eliminated NIR-excess sources from the sample using 
the following algorithm. First, we fit a vector to the data (a straight line in the color-color 
diagram) using all the stars with photometric errors of less than 0.01 mag. This vector will 
have a flatter value of Ej-h/Eh-k than the true reddening vector since stars with excess 
are located below it. Then, we remove all the stars having NIR-excess with respect to this 
vector and fit a new vector to the remaining stars. We repeat this process until the value of 
Ej_ H /E H _ K converges. Figure \5\ shows the resulting reddening vectors for the four regions. 
The values of Ej-h/ Eh-k range between 1.71 and 1. 86 and are in agreement with previou s 
values from the literature (Ej^h/Eh-k = 1-7 to 2.1. 1Kenyon et al.lll998l ; Ijiang et al.ll2002l ). 
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3.2. Color-magnitude diagrams 

Color-magnitude diagrams for all star forming regions are presented in Figure [61 NIR- 
excess sources are marked with asterisks. We used the color-magnitude diagram to divide 
the population of stars in two groups: massive stars (stars with colors corresponding to 
spectral type earlier than B2) and intermediate-mass stars (with colors corresponding to 
spectral types between B2 and A5). The magnitude of an unreddened A5 star located at 
the same distance than G329.337 (the furthest region) is K=16 mag, significantly brighter 
than our completeness limit (see Table HJ) indicating that we are complete in detections for 
both groups. The total number of stars earlier than A5 detected in each region is listed in 
Table EJ 

Several sources were detected only in the H and K bands (see Table [3]). These J- 
dropout sources are likely the most embedded YSOs. However, it is difficult to distinguish 
YSOs from main sequence stars among the J-dropout sources because in this case the color- 
color diagnostic used to find NIR-excess sources cannot be used. Instead, we selected J- 
dropout sources that are brighter than a B2 type star. Since the K-band is more sensitive 
to IR-excess, YSOs will be preferentially brighter in this band. The location of the selected 
J-dropout sources in the H-K vs. K color-magnitude diagram are indicated by the squares 
in Figure EJ Their spatial distribution is discussed in § 13.41 



3.3. Surface density of NIR-excess sources 

Having identified the sources with NIR-excess, we can study their spatial distribution. 
This was done by dividing the mosaics using a regular grid and calculating the surface density 
of stars at each point of the grid: 



where is the distance to the N = 10 nearest neighbor (NN). Figure [7] shows the surface 
density maps of NIR-excess sources for all the regions. The maps were created using a 3 
arcsecond grid and smoothed applying a convolution with a 3.0 arcsecond Gaussian kernel. 
The lower contour in the maps correspond to 3<r over the mean value. The peak surface 
densities are 120, 40, 75 and 30 stars per pc 2 for G324.201, G328.307, G329.337 and G330.949, 
respectively. 

We use the surface density maps to find groups of NIR-excess sources. We define a group 
as a concentration of 10 or more NIR-excess sources with a surface density of stars higher 
than 3a over the mean. Using this definition, we found five groups: two towards G324.201 
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and G328.307 (labeled A and B) and one towards G329.337. Spurious group detections near 
the edges of the field in each region were ignored due to large photometric errors. Table [5] 
gives group parameters. Columns 2 and 3 list the coordinates corresponding to the peak 
surface density. Figure [7] shows also the spatial distribution of massive and intermediate- 
mass sources with NIR-excess (in red and blue respectively). Groups G324.201A, G328.307A 
and G329.337 contain a large number of massive NIR-excess sources while G324.201B and 
G328.307B are mainly composed by intermediate and low-mass stars. The peak position of 
groups G324.201A, G328.307A and G329.337 is coincident, within the errors, with the peak 
position of the massive dense cores detected at 1.2 mm in Paper I. Groups G324.201B and 
G328.307B are found projected towards the massive and dense cores and are less prominent 
than their A companions. In what follows, we consider groups A and B as substructures 
of a single cluster, centered at the peak position of group A. The total number of cluster 
members, estimated by counting the NIR-excess sources inside a radius of 1.5 pc from the 
cluster center, is given in Table El This radius corresponds roughly to the radius of the 10% 
level of the peak 1.2 mm dust emission from the massive and dense cores (Paper I). 

Since the NIR excess emission is likely to arise from disk-like structures, we assume that 
the NIR-excess sources are young and hereafter we will refer to them as YSOs. The spatial 
distribution of massive and intermediate-mass YSOs is discussed in § 14.11 The non-detection 
of YSO groups in G330.949 is discussed in § 14.51 

3.4. Spatial distribution of J-dropout sources 

As mentioned in § 13.21 the J-dropout sources selected (brighter than a B2 star) likely 
correspond to embedded YSOs. If massive J-dropout sources correspond to either field stars 
or background AGB galaxies, we expect them to be uniformly distributed over the field. On 
the other hand, if they are associated with young clusters of NIR-excess sources discussed 
in § 13.31 we expect them to have a similar spatial distribution. 

The surface density maps of massive J-dropout sources are shown in Figure El They were 
created using the same technique as for the NIR-excess surface density maps. Figure [8] shows 
that far from being distributed homogeneously, the sources are concentrated at the location 
of clusters G324.201, G328.307 and G329.337, suggesting that those J-dropout sources are 
cluster members. These are probably very embedded sources and good protostar candidates. 
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3.5. Extinction maps 



We compare the distribution of young stars with the distribution of dust. We created 
extinction maps in each region using the mean value of the K-band extinction (Ak) of the 
nearest N = 10 stars over a uniform 3.0 arcseconds grid. The Ak values were estimated by 
derredening the sources in the color-color diagram to their main sequence colors. Sources with 
NIR-excess were derredened to the CTTS loci (see Figure HJ). We convolved the resulting map 
with a 3.0 arcsecond Gaussian kernel. In order to eliminate the contribution of foreground 
stars to the extinction we used stars with Ak > 1.0 in G324.201 and G329.337 and with 
Ak > 0.8 in G328.307 and G330.949. This corresponds to the interstellar extinction at the 
distance to each region a ssuming an interstellar extinction of Ak = 0.15 magnitudes per kpc 
( llndebetouw et al.l 120051 ). It is important to note that due to the probable non detection of 
background stars located behind densest parts of the molecular clouds, the extinction values 
estimated represent a lower limit to the true extinction towards those regions. 

Figure M shows the K-band extinction map in each of the four regions. Groups G324.201 A 
and G328.307A coincide with extinction peaks. The extinction peak in G330.949 is located 
at less than 1 arcsecond from the peak 1.2 mm emission (indicated by the arrow in Fig. [9]). 
The mean values of Ak are close to 1.5 in all regions and the peak Ak values are 2.2, 2.7, 2,2 
and 2.6 for G324.201, G328.307, G329.337 and G330.949 respectively. These values are sim- 
i lar to the extinction i n other sites of hi gh mass star format ion like the complex S254-S258 
dChavarria eiTallhoOsh and DR21/W75 fcumar et al1l2007h . 



4. Discussion 

4.1. Radial distribution of NIR-excess sources inside the massive and dense 

cores 

In § 13.11 and 13.21 we used the color-color diagram to identify NIR-excess sources and 
the color-magnitude diagram to divide stars into two groups depending on their estimated 
spectral type: massive and intermediate-mass stars. In this section, we investigate the spatial 
distribution of the NIR-excess sources from both groups within the massive and dense cores. 
We define the radial surface density of embedded sources as: 

Ni 



where iVj is the number of NIR-excess sources located between a radius of r»_i and r, from 
the center of the cluster. The step in radius used is 0.2 pc, about a third of the radius of 
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the massive dense cores. Figure [10] shows £, histograms for massive and intermediate-mass 
YSO. The error bars in Figure [10] correspond to the square root of the number of YSO 
in each bin. The histograms show a clear difference in the radial distribution of massive 
and intermediate-mass YSO. In clusters G328.307 and G329.337, massive YSO are highly 
concentrated within the central region whereas intermediate-mass YSO are more spread. 
The mean distance to the cluster center for massive YSO in G328.307, (r^), is 0.09 pc and 
for G329.337 is 0.15 pc, while for intermediate-mass YSO the mean distances (r : ) are 0.60 
and 0.87 pc, respectively. 

Cluster G324.201 shows a concentration of both massive and intermediate-mass YSO 
in the core. To quantify how different are both distributions, we fitted the surface density 
S with a radial power-law of the form: 




where r is a constant. The values of a in G324.201 for massive and intermediate-mass YSO 
are 1.79 ± 0.06 and 0.77 ± 0.09, respectively, suggesting that massive IR-excess sources are 
more concentrated than intermediate-mass YSO. 



4.2. Dynamical effects and cluster ages 



In the previous section we concluded that there are signatures of mass segregation in the 
clusters G324.201, G328.307 and G329.337. In order to determine whether the segregation 
is primordial or dynamic (see § [TJ, we will compare the age of the cluster members with 
their relaxation tim e. One of the most used expressions to estimat e the relaxation time of 
clusters is given by (ISpitzer fc Hartlll97ll ; iBinney fc Tremaindll987l ); 

^^(^^(^('Sfyr, (4) 



ln(0.4iV) V.10 5 M G 



ni 



where N is the number of cluster members, M is the total mass, m is the mean mass 
of cluster members and r H is the half-mass radius (approximately 0.5 pc for all clusters). 
We calculated the number of cluster members by counting the YSOs inside a radius of 
1.5 pc from the cluster center. The number of undetecte d cluster mem b ers w as estimated 
using a synthetic cluster created with the algorithm from iMuench et al.l (120001 ) . We used a 
Trapezium IMF and age of 0.5-1.0 Myr as inputs. Based on the number of observed YSO, we 
used the synthetic cluster to estimate the percentage of cluster members with a magnitude 
fainter than the K-band completeness limit (see Tabled]). The corrected number of cluster 
members is then estimated by correcting the observed number of NIR-excess sources by the 
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percentage of undetected members. Having the corrected number of cluster members, their 
mass is estimated assuming a mean mass of m — 0.5 M per star. Using equation (j3J), we 
found relaxation times t re i ax between 5 and 8 Myr for cluster members in the studied dense 
cores. 

In order to estimate the age of the clusters, we can use for instance the fact that they 
are still associated with their parental molecular cloud, implying that they are younger than 
5 Myr (ILeisawitz et al.lll989l ). A more accurate estimate can be obtained using theoretical 
isochrones in the color-magnitude diagram. We consider NIR-excess sources within 1 pc 
from the center of the clusters and compare their positio n in the color-mag nitude diagram 
with the 1 Myr pre-main sequence (PMS) isochrone from ISiess et al.l (120001 ) . To reproduce 
the cluster environment we added an extinction of = 2.0 to the isochrone. Figure [TT] 
illustrates this method for G328.307, showing that cluster members have ages of less or of 
the order of 1 Myr. Using the same method, the estimated ages of cluster members for 
G324.201 and G329.337 are also of the order of 1 Myr. This suggests that the mentioned 
stellar clusters are not relaxed yet and that equipartition of energy is not the origin of the 
observed mass segregation. 

Since the dense cores are very young and mainly gaseous, with only a small percent- 
age (~10%) of the total mass in the form of stars (see Table [6]), dynamical friction (also 
called gravitational drag) produced by the gaseous background onto the stars may play an 
important role in the orbital evolution of stars. Dynamical friction is defined as momen- 
tum loss by a massive perturber (in this case a star) due to the interaction with it s own 
gravitationally induced gaseous wake. It was first studied by IChandrasekharl (119431 ) for a 
background medium composed only by sta rs, but it can be extended to the evolution of a 
mass i ve perturber in a g aseous background (IOstrikerlll999l ; ISanchez-Salcedo fc Brandenburg 



200ll ; iKim fc Kiml 120071 ). In this case, we find that for a system of density p and Mach 



number . M, a massive p erturber of mass M star and velocity V star feels a gravitational drag 
given by (IOstrikerlll999l ): 



DF 



-47rp 



GM. 



star 



star 



xf(M) 



(5) 



where f (M.) is a dimensionless funtion than depe nds, in additio n to the Mach number, on 
the maximum and minimum impact parameters (jOstrikerl Il999l ). This formulae has been 
numerically tested, finding surprinsingly good agreement with numerical resu lts when the 
maxi mum impact parameter equals to the diameter of the perturber's orbit (IKim fc Kim 
2007h . 



Gravitational drag on a gaseous background has been applied to study the evolution 
and migration of stars in dense gaseous star forming clouds using numerical smooth particle 
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hydrodynamics simulations (jEscala et al.ll2003l 12004 , in § 5.2). The simulations have initial 
conditions similar to the ones found in the dense cores studied: the cloud density profile is 
described by an isothermal sphere (p(r) oc r _L8 ), with perturber having masses of 1 % of 
the total gas mass (M star ~ 20 M « 0.01 M gas ). The characteristic ti mescale found in th ese 



simulations for migration of the stars by a factor of ~2-3 is given by (jEscala et al.ll2004 , eq. 
18): 



T 



M 



10 M 



-1/2 



R 



0.01 pc 



3/2 



3.6 x 10 4 yr, 



(6) 



where M gas is the total mass of the gas and R is the initial distance of the stars to the center 
of the cloud. 

Using equation (6) for a typical cloud mass of 4 x 10 3 M and initial distance of 0.4 pc, 
the migration timescale is r = 4.5 x 10 5 yr. This is comparable to the estimated ages of 
the mentioned clusters (~ 1 Myr). Furthermore, the dependence of the gravitational drag 
force on the perturber's mass (equation 5), gives migration timescales proportional to the 
inverse of the perturber's mass. This implies that intermediate mass stars have migration 
timescales typically a factor ~10 longer than those of high mass stars and hence longer than 
the estimated age of the clusters (~ 10 6 yr). Therefore, gravitational drag is considerably 
less effective for intermediate mass stars compared with the most massive ones. This suggest 
that in the mentioned stellar clusters only the massive stars were able to efficiently migrate 
towards the center due to gravitational drag, being a possible dynamical explanation for the 
origin of the observed mass segregation. 



4.3. Sample incompleteness and clusters IMF 



The analysis performed in § 14.11 sh ows evidence of mass segregation in three massive and 
dense cores. In a recent paper, lAscenso et al.l (120091 ) studied the effects of binning and data 
incompleteness on the most common diagnostics used to find mass segregation. They found 
that incompleteness due to crowding in the center of non-segregated clusters will mimic the 
observed evidence of mass segregation. This is because low-mass stars located near massive 
stars are blended from detection by the higher luminosity of their massive companion and 
remain undetected, producing a fictitious lack of low-mass stars in the center of the cluster. 

Since th e clusters stud i ed in this paper have different properties than the synthetic 
clusters from lAscenso et al.l (120091 ). we can't use their results in order to correct the incom- 
pleteness of our data. Thus, we investigated how the different sources of incompleteness - 
dust extinction, number of cluster members and crowding - may affect our results. The 
effects of extinction are difficult to evaluate since extinction is unique for every star. The 
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estimated number of cluster members (see §3.31) corresponds to a lower limit s i nce w e are not 
taking into account possible members without NIR-excess. IChavarria et al.l (120081 ) showed 
that the percentage of non-excess members in young clusters range roughly between 20- 
80% for clusters with ages between 0.5-2 Myr. Since the clusters studied here are probably 
younger than 1 Myr, we expect that the percentage of non-excess members is roughly 50% 
and therefore we are identifying approximately half of the cluster members. 



Ascenso et al.l (120091 ) showed that crowding will have an important effect near the center 



of the cluster. The identification of two or more stars as separate entities near the central 
region will be more difficult due to the small distance between cluster members. Two stars 
will be identified as separate entities if their projected distance is about two times the angular 
resolution. For our data, this translates roughly in a separation of 0.02 pc at the distances 
of the clusters (see § 12.31) . Assuming, for simplicity, that the stars near the cluster center are 
distributed homogeneously, then we will be able to identify separate entities if the density 
of stars is less than 2.5 x 10 3 stars per pc 3 . This value is higher than the density of massive 
plus intermediate-mass stars in the central parsec of the dense cores (~ 150 stars per pc 3 ) 
even if we correct by the number of sources with no IR-excess, indicating that crowding is 
not an important source of incompleteness. 

A possible way to test the evidence of mass segregation in our data is to calculate, for 
a given IMF, the number of intermediate-mass YSOs expected from the high-mass YSOs 
detected and compare it with the observed number. If both numbers agree, the data is 
complete and mass segregation is real. We define the IMF by: 



dN 



d(logm) 



Am 1 



(7) 



where N is the number of stars, m is the mass, T is the IMF slope and A is a constant. 
From the number of high-mass NIR-excess sources, we estimate the number of i ntermediate- 
mass stars using two different IMF slopes: Salpeter (r = —1.35, IScalol Il986l ) and Arches 
(r = —0.9, IStolte et al.ll2005l ). The mass range used for the high-mass stars is 30 to 8 
M and for the intermediate- mass is 8 to 2.1 M . Montecarlo simulations were used to 
estimate the probability of drawing Nim intermediate-mass stars given the number Nm of 
massive stars from a mother distribution function described by a power-law with exponent 
T = —1.35 and —0.90. Table [7] shows that the distribution of mass in G324.201 is better 
described by an Arches IMF, in G328.307 the distribution of mass is well described by 
both Salpeter and Arches IMF. This suggests that mass segregation in these clusters is real 
assuming a top-heavy IMF. 



The hypothesis of a top-heavy IMF is also suggested by the comparison of the gas mass 
with the mass in stars for different IMFs in dense cores. The mass of gas in dense cores 
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was estimated in Paper I using two different methods: the CS emission line (M cs ) and the 
1.2 mm dust emission (M du8t ). Both masses are in excellent agreement for all the dense cores. 
In addition, the total mass of the cores (gas plus stars, M virial ) was also estimated, obtaining 
similar values to the mass of gas in cores. This means that the mass of stars is negligible 
compared to the mass of gas for most of the cores, which is also suggested by our estimation 
of the mass of stars (see Table [6]). If we use the number of cluster members with masses of 
more than 8 M to estimate the total mass of stars for different IMF slopes, we find that 
for a Salpeter IMF, the mass of stars is about half of the mass of gas and non negligible (see 
Table [Sj). This contradiction suggest that there is a lack of low-mass stars in the studied 
clusters and supports the idea of a top-heavy IMF. 



4.4. Different main sequences in color-magnitude diagram 



In this section we investigate the origin of the multiple main sequences seen in the color- 
magnitude diagrams of G328.307, G329.337 and G330.949 (see Figure [6}. All the studied 
sources are located near the tangent of the Norma spiral arm. along those lines of sight 
there are two other spiral arms at closer distances: Carina and Scutum-Crux. Molecular 
clouds that belong to those spiral arms are likely to cover the field of view adding layers of 
dust. They will increase the extinction of background stars which will appear as reddened 
main sequences in the color-magnitude diagram. Figure [T2] shows three main sequences at 
different extinction in the color-magnitude diagram of G330.949 (NIR-excess sources were 
plotted as points to distinguish the different main sequences). 

We correlate the column density derived from the extinction of the different main se- 
quences in the color-magnitude diagrams with the column d ensity derived from mo lecular 
material in the line of sight to the clusters. The 12 CO survey of iBronfman et al.l (119891 ) shows 
that the clusters have a Vlsr of the order of -120 to -80 km s^ 1 , but it also shows that there 
is molecular material with Vlsr > ~ 80 km s _1 along the line of sight to all the regions. In 
the case of G330.949, there are two molecular components; one with Vlsr ~ — 65 km s _1 
and other with ~ —45 km s _1 . We esti mated the H2 co lumn density from the extinction in 
the K-band using the relation given by iDickmanl (119781 ): 



N(H 2 ) 
A v 



1.25 x 10 21 cm 2 mag~ 



and the relation between visual and K-band extinction from ICardelli et al.l (119891 ): 

^ = 0.114. (9) 

Ay 
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By combining equations (jHJ) and © we obtain: 

N ^ ~ 1.1 x 10 22 cm" 2 mag' 1 . (10) 



.4 



A' 



We used the As: values of the main sequences in the color-magnitude diagram (see Fig- 
ure [T2j) and equation (fit)]) to estimate the H2 column density. The values estimated are 
0.58xl0 22 cm" 2 for A K = 0.55 and 0.48xl0 22 cm" 2 for A K = 0.45. On the other hand, 
H9 column dens i ty val ues derived from the molecular component were taken directly from 



Bronfman et al.l (119891 ). being in agreement by a factor of 2 with the column densities de- 
rived from Ak- A comparison of the H2 column density values derived from both methods 
is shown in Table [9] for G328.307, G329.337 and G330.949. We conclude that the reddened 
main sequences in the color-magnitude diagrams are probably due to the presence of layers 
of dust along the line of sight. 



4.5. The stellar population within the G330.949 molecular core 



As discussed in § 13.31 the G330.949 molecular core is the only core in which we did not 
detect a group of NIR-excess sources, even though this region is the most luminous at far-IR 
wavelengths (see Tabled]). To investigate the main sources of luminosity in this region we 
examined individual objects within the contours of the 1.2 mm emission. We found seven 
objects, labeled in Figure [T31 with characteristics commonly associated to star formation: 
NIR-excess, extended emission and high extinction. Objects 2 and 3 are the most extinted 
stars in the region as can be seen in the color-color diagram (Fig. H]). Objects 1 and 4 are 
J-dropouts and object 5 has NIR-excess. In addition, objects 1 and 5 are associated with 
extended emission. Most of these objects are concentrated near the peak position of the 
1.2 mm emission suggesting that star formation is taking place near the center of G330.949. 
The position of objects 1, 2 and 3 in the (K vs. H-K) color-magnitude diagram suggests that 
they correspond to highly extincted 06 stars. If so, their total luminosity is ~ 7 x 10 5 L , 
indicating that they are responsible for most o f the luminosity of th e region. Objects 1, 2-3 



and 6-7 correspond to objects 1, 2 and 3 from lOsterloh et al.l (ll997l ). They identified these 



objects as highly reddened sources in agreement with our results. 

In addition, radio continuum observations with high (~ 1") angular resolution towards 
IRAS 16060-5146 showed the presence of a co mpact HII region w ith an irregular morphology 



and a peak flux density of 280 mJy at 8.6 GHz (IWalsh et al.lll998l ). It lies at the peak position 



of the dust core and is not associated with any of the NIR sources in our images, suggesting 
that the exciting source of the H11 region is extremely embedded and undetectable at NIR 
wavelengths. 
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5. Conclusions 

We used NIR data from the VLT-ISAAC instrument to identify the young population of 
stars in four regions of massive star formation: G324.201, G328.307, G329.337 and G330.949. 
These regions are located in the Norma spiral arm at distances between 5.4 and 7.3 kpc from 
the Sun. We used color criteria to identify sources with NIR-excess and found 236, 120, 186 
and 228 NIR-excess sources in G324.201, G328.307, G329.337 and G330.949 respectively. 

We identified clusters of YSOs in surface density maps of NIR-excess sources. We found 
new clusters towards G324.201, G328.307 and G329.337. The spatial distribution of sources 
detected in H and K, but not in the J-band (J-dropout), was investigated. We concluded that 
J-dropout sources with colors of massive stars are likely to correspond to cluster members. 

We analyzed the spatial distribution of massive and intermediate-mass cluster members 
and found evidence of mass segregation in each of the clusters. The effects of data incom- 
pleteness on the detected mass segregation were assessed and concluded that they are not 
important. 

The age of the clusters was compared with their relaxation time and migration timescales 
due to gravitational drag, in order to determine if mass segregation is primordial or dynam- 
ical. We found that the migration timescales for massive stars is comparable with the age 
of the cluster members. This evidence suggest that gravitational drag from the gaseous 
component may be a possible cause for mass segregation in clusters G324.201, G328.307 and 
G329.337. 

Several main sequences in the color-magnitude diagram of regions G328.307, G329.337 
and G330.949 were identified. We computed the H 2 column density derived from 12 CO for 
molecular components located in Scutum-Crux arm and compared those values with the 
column density derived from for the different main sequences. The H2 column densities 
derived from both methods are in reasonable agreement, indicating that the reddening is 
mainly due to molecular clouds located along the line of sight. 

This work was supported by CONICYT-CHILE through projects FONDAP 15010003 
and BASAL PFB-06. 

Facilities: VLT(ISAAC). 
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Fig. 1. — Errors in the K-band photometry for G324.201. Upper panel: error as a function 
of K-band magnitude. Lower panel: error as a function of the distance to the center of the 
mosaic. The segmented line indicates the maximum error allowed. The error is independent 
of the distance to the center of the mosaic. 




Fig. 2. — JHK-color image of the regions observed (J-blue, H-green, K-red). The white bar 
represents 1 parsec. The arrow indicates the position of the center of the clusters. 
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Fig. 3. — Apparent magnitude histogram for detections with photometric errors < 0.2 magni- 
tudes in all bands for G328.307. The vertical lines on the top indicate the 90% completeness 
limits in the J, H and K bands: 20.7, 18.5 and 17.5 mag respectively. 
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Fig. 4. — JHK color-color diagrams for the four regions. Gray lines correspond to the 
position of the main sequence (light gray) and giant stars. Dashed lines enclose reddened 
main sequence and giant stars. The dashed-dotted lines are displaced by the A dispersion of 
the sources located above the giant branch. S ources with NIR-ex cess are marked as asterisks. 
The dotted lines indicate the CTTS loci from lMeyer et al.l (119971 ). The percentage of sources 
located above the giant branch is 3.5, 10, 6 and 2.5% for G324.201, G328.307, G329.337 and 
G330.949 respectively. Gray squares in G330.949 correspond to individual sources from 
Figure [Lj (see §03]). 
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Fig. 5. — JHK color-color diagram used to determine the extinction law in each region. We 
used stars with photometric errors smaller than 0.01 magnitudes in all bands. Points marked 
as "x" correspond to stars eliminated by the algorithm explained in § 13.1.11 The solid lines 
indicate the extinction law found in each region. 
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Fig. 6. — Color-magnitude diagrams for the regions observed. The gray lines show the 
position of the main sequence stars (jKoornneeil Il983l ) at the adopted distances. Dashed 
lines correspond to the reddening vector for BOV, B2V and A5V stars. Sources with NIR- 
excess are marked as asterisks. Selected J-dropout sources are marked as empty squares. 
The diagrams show that the color of the stars is shifted from the main sequence due to 
the interstellar extinction. It is possible to identify several main sequences with different 
foreground extinction in each region. Black squares in G330.949 correspond to the individual 
sources from Figure [13] (see § 14.51) . 
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Fig. 7. — Surface density maps of NIR-excess sources. White contours begin at 3cr over the 
median surface density (22, 11, 13 and 12 stars per pc 2 for G324.201, G328.307, G329.337 
and G330.949 respectively) and increase by 3a. Massive and intermediate-mass stars with 
NIR-excess are marked with red and blue dots, respectively. Black contour levels correspond 
to 10, 30, 50, 70 and 90 % of the peak dust emission detected with SIMBA at 1.2 mm 
( iGaray et al.l 120091 ) . The straight line indicates a length of 1 pc. 
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Fig. 8. — Surface density maps of massive J-dropout sources with colors corresponding to 
spectral types earlier than a B2V star (shown as squares in Figure[6]). Black contours begin at 
3a over the mean surface density (9.3, 3.3, 15 and 2.25 stars per pc 2 for G324.201, G328.307, 
G329.337 and G330.949 respectively) and increase by 3cr. White contours correspond to the 
surface density of NIR-excess sources from Figure [71 The straight lines represent a distance 
of 1 pc. There is a clear spatial correspondence between massive J-dropout and NIR-excess 
sources in G324.201, G328.307 and G329.3 37. The white plus-s ign in G330.949 indicates the 
position of the 1.2 mm emission peak from iGaray et al.l (120091 ) . 
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Fig. 9. — K-band extinction maps for the observed regions. Contours begin at 1.8 A# and 
increase by 0.2. The black lines represent a distance of 1 pc. The center of the clusters 
are indicated by th e arrows. In G330.949, the arrow indicates the position of the 1.2 mm 
emission peak from I Gar ay et al.l (120091 ) . 
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Fig. 10. — Radial surface density of NIR-excess massive (gray) and intermediate-mass (line) 
sources in clusters G324.201, G328.307 and G329.337. 
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Fig. 11. — Color-magnitude diagram for G328.307. Stars correspond to NIR-excess sources 
within 1 pc from the center of the cluster G328.307. The dot-dashed line indicates the 1 Myr 
pre-main sequence track from lsiess et al.l ( 2000 ). We added an extinction of Ak = 2.0 to the 
PMS model to resemble the physical conditions of the cluster. The dashed line corresponds 
to the reddening vector for a K7 PMS type star. The arrow corresponds to an excess of 0.5 
mag in the H — K color, which is the mean excess of cluster members in G328.307. 
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Fig. 12. — Color-magnitude diagram for G330.949. We only included sources detected in J, 
H and K bands. Arrows indicate the extinction of the reddened main sequences at of 
0.55, 1.0 and 1.9. 
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Fig. 13. — K-band image of G330.949. Seven sources of interest are shown in the figure (see 
explanation in § 14.51) . Contours correspond to 10, 30, 50, 70 and 90% of the 1.2 mm dust 
emission (from Figure [7]). The black line represents 1 pc. 
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Table 1. Source list 



Source 


IRAS id. 


Ion 


lat 


a(J2000) 


<5(J2000) 


d 


Vlsr 8, 


LlRAS 






(deg) 


(deg) 


hh mm ss 


dd mm ss 


(kpc) 


(km s- 1 ) 


(Ls) 


G324.201 


15290-5546 


324.201 


0.119 


15 32 52 


-55 56 19 


6.9 b 


-88.0 


7.8 x 10 5 


G328.307 


15502-5302 


328.307 


0.432 


15 54 05 


-53 11 37 


5.6 C 


-91.7 


1.1 x 10 6 


G329.337 


15567-5236 


329.337 


0.147 


16 00 32 


-52 44 52 


7.3 b 


-107.1 


1.3 x 10 6 


G330.949 


16060-5146 


330.949 


-0.174 


16 09 52 


-51 54 38 


5.4 C 


-89.7 


1.0 x 10 6 



a From iBronfman et alJ (|1996l ). 

b Tangent distance to Norma arm. 



"Near side distance to Norma arm (jBronfman et al.ll2000l ). 
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Table 2. Integration Times 



Source J s a H a K s a 



G324.201 1200 2400 3600 

G328.307 1200 2400 3600 

G329.337 300 300 300 

G330.949 300 300 300 



a Time in seconds. 
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Table 3. Number of detected stars a 



Source J s H K s J-dropout b J S HK S M c IM d 



G324.201 4809 5721 8347 

G328.307 4103 6036 7057 

G329.337 4431 6357 7927 

G330.949 3516 5512 6887 



1081 4057 99 1449 

1349 3779 64 1002 

1424 3876 82 1411 

1580 2894 49 773 



a Detections with error < 0.2 mag. 

b Sources detected in H and K bands but no detected in J band. 
c Massive: stars with spectral type earlier than B2. 
d Intermediate-mass: stars with spectral type between B2 and A5 
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Table 4. 90% Completeness 



Source 


Js 


H 


K s 




(mag) 


(mag) 


(mag) 


G324.201 


20.2 


17.5 


17.2 


G328.307 


20.7 


18.5 


17.5 


G329.337 


19.5 


17.5 


16.5 


G330.949 


19.3 


16.8 


16.7 
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Table 5. Group properties 



Name a(J2000) <5(J2000) N m a N M b N M C 
hh mm ss dd mm ss 



G324.201A 


15 


32 


52.8 


-55 


56 


15.7 


28 


12 


9 


G324.201B 


15 


32 


49.3 


-55 


55 


58.2 


15 


1 


5 


G328.307A 


15 


54 


06.5 


-53 


11 


39.9 


20 


7 


6 


G328.307B 


15 


54 


06.6 


-53 


12 


02.0 


10 





5 


G329.337 


16 


00 


33.1 


-52 


44 


45.6 


26 


12 


5 



a Number of NIR-excess sources. 
b Number of massive NIR-excess sources. 
c Number of intermediate-mass NIR-excess sources. 
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Table 6. Cluster parameters 



Cluster 




N b 

1,1 cor 


Mstars 


M dust c 








[M ] 


[M ] 


G324.201 


130 


890 


445 


5.0xl0 3 


G328.307 


60 


345 


173 


5.6xl0 3 


G329.337 


62 


1240 


620 


5.8xl0 3 



a Number of NIR-excess sources inside a ra- 
dius of 1.5 pc from the cluster center. 

b Corrected number of members (see § 14.21) 

c From Paper I 
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Table 7. Expected number of intermediate-mass members in dense cores 



Cluster 


N M a 


NiM b 


NIM Salp eter C 


r> d 

P-1.35 


NIM Arches e 


P-0.90 


G324.201 


20 


55 


122 


0.03 


67 


0.49 


G328.307 


7 


30 


43 


0.43 


23 


0.51 


G329.337 


18 


19 


110 


0.00 


60 


0.02 



a Number of massive IR-excess sources. 
b Number of intermediate-mass IR-excess sources. 
c Expected intermediate-mass stars for T = —1.35. 
Probability of drawing N IM given N M for a T = —1.35 IMF. 
c Expected intermediate-mass stars for T = —0.90. 
f Probability of drawing N IM given N M for a T = -0.90 IMF. 
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Table 8. Mass of stars and gas in dense cores 



Cluster 


N a 


Msalpeter b 


M Arches° 


M cs d 


M dust d 


■^-virial 






[M ] 


[M ] 


[M ] 


[M ] 


[M ] 


G324.201 


43 


3000 


1700 


4.3xl0 3 


5.0xl0 3 


3.7xl0 3 


G328.307 


43 


3000 


1700 


4.6xl0 3 


5.6xl0 3 


6.8xl0 3 


G329.337 


60 


4100 


2300 


6.0xl0 3 


5.8xl0 3 


6.0xl0 3 


G330.949 








l.OxlO 4 


l.lxlO 4 


7.9xl0 3 



a Number of massive IR-excess plus J-dropout sources inside a 1.5 pc 
radius. 

b Mass of stars between 0.6 and 30 M for T = -1.35. 
c Mass of stars between 0.6 and 30 M for T = —0.90. 
d From Paper I. 
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Table 9. H 2 column density derived from 12 CO and k K . 



Source 


i2 CO a 




i2 CO b 


A b 
A K 




10 22 [cm- 2 ] 


10 22 [cm- 2 ] 


10 22 [cm- 2 ] 


10 22 [cm- 2 ] 


G328.307 


0.78 


0.53 


0.30 


0.65 


G329.337 


0.73 


0.85 


0.50 


0.65 


G330.949 


1.05 


0.58 


0.57 


0.48 



a Component with Vlsr between -30 and -57 km s 
b Component with Vlsr between -57 and -80 km s 



